Demonstrating various quantum effects with two entangled laser beams 
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We report on the preparation of entangled two mode squeezed states of yet unseen quality Based 
on a measurement of the covariance matrix we found a violation of the Reid and Drummond EPR- 
criterion at a value of only 0.36 dz 0.03 compared to the threshold of 1. Furthermore, quantum state 
tomography was used to extract a single photon Fock state solely based on homodyne detection, 
demonstrating the strong quantum features of this pair of laser-beams. The probability for a single 
photon in this ensemble measurement exceeded 2/3. 



Introduction. Entanglement plays an ubiquitous role 
in the field of quantum physics [IJ. Light fields are 
w^idely used as the carrier for entangled states as they 
are readily accessible and of high quality in the experi- 
ments. Traditionally, quantum opticians are in two com- 
plementary camps. One is dealing with discrete vari- 
ables using individual photons, the other with contin- 
uous variables using bright laser beams. Both camps 
have collected a rich set of tools and techniques in 
their individual toolboxes. On the one hand, the dis- 
crete camp focuses on the preparation and counting 
of photons. The corresponding theory describing the 
quantum states and the corresponding measurements is 
based on finite dimensional Hilbert spaces. Experimen- 
tal progress has included various sources for single (and 
multiple) photons[2. 3|, sources for entangled pairs and 
larger groups of entangled photons IK O. Numerous Bell 
tests, quantum teleportation, different aspects of quan- 
tum communication were demonstrated! 6 -8 1. On the 
other hand, the continuous camp focuses on the field 
aspects of quantum light and uses mainly phase sen- 
sitive homodyne detection. In order to accomplish an 
equivalently simple theory the theorists had to stick to 
Gaussian fluctuations and operations. This is well justi- 
fied for the majority of present experiments which have 
access to moderately strong nonlinearities. There are 
squeezed states, Einstein-Fodolsky -Rosen (EPR) like en- 
tangled states, and an equal variety of quantum commu- 
nication topics! 9 -12 J. A common aim here is to increase 
the strength and reliability of the entanglement in order 
to improve the quality of the applications. 

In recent years hybrid experiments have overcome the 
separation of discrete and continuous variable quantum 
optics e.g. !'3l [T3HT6l . These combine single photon and 
homodyne detection. The cZ/cfc-detector is used to gen- 
erate trigger pulses displaying a non-Gaussian opera- 
tion. These trigger pulses are used to gate the homodyne 




Figure 1: Wigner function of a single photon state recon- 
structed from a pair of entangled laser-beams. We used quan- 
tum state tomography in order to characterise the output state. 
The reconstruction shown here is based on the inverse radon 
transform without any correction for detection efficiency or 
detector dark noise. 



measurement according to the temporal shape of the de- 
tected photon, which then produces states of light with 
impressively strong quantum features. Based on photon 
pair sources this method was used to generate e.g. single 
(few) photon states, single photon added thermal states 
and single(few) photon subtracted squeezed states! 3, 
[T3tiT7l . The latter exhibit a high fidelity with small 
Schrddinger cat states. The goal is to show these effects as 
clearly as possible. However, present single photon de- 
tectors suffer either from poor mode (spatial, frequency, 
temporal) selectivity, limited detection efficiency, notice- 
able dark noise, lacking photon number discrimination 
or a combination of these. 



The work we present here entirely relies on contin- 
uous techniques, i.e. we use a bright pair of quadra- 
ture entangled laser beams and balanced homodyne de- 
tection only. A number of technical improvements to 
our experimental apparatus lead to a previously unseen 
quality of EPR entangled states as described by the well 
known criterion for bipartite continuous variable entan- 
glement |[T8ti2Ql . For the combined conditional vari- 
ance we found a value as low as 0.36. Additionally, 
we demonstrate the proposal of Ralph, Huntington and 
Symul [21 J reconstructing a single photon Wigner func- 
tion by post processing our experimental data, which 
is analogues to the hybrid experiments without actu- 
ally detecting single photons. This method is similar 
to the generation of heralded single photons as based 
on spontaneous parametric downconversion, and was 
also used for the reconstruction of optical Schrodinger 
Kitten states [22J. However, instead of using a single 
photon detector on one of the output modes we used 
homodyne detection in both and exploited the relation 
n + 2 = i + ^^)/ with n representing the photon 
number, X and P for the quadrature operators. This idea 
was suggested in ||23| and applied in [24J for the first 
time. Of course, we cannot measure X and P simulta- 
neously, but rather only sequentially. However, relying 
on our ability to prepare the same state over and over 
again, in the ensemble average we can use the relation 
in order to mimic the real photon detection. The recon- 
structed probability for a single photon for the experi- 
mental data of Figure |4] exceeds 67 % and can be seen as 
an overall figure of merit for the efficiency. Note that we 
did not correct our data for any experimental imperfec- 
tions, such as quantum efficiencies or dark noise unless 
indicated otherwise. 

Preparation of Entanglement. A scheme to generate a 
pair of entangled beams is shown in Fig. |2] The main 
building block is a pair of optical parametric amplifiers 
(OPA) operating in a de-amplification regime such that a 
pair of amplitude-squeezed states of light is produced. 
A dual wavelength laser at 1064 nm and 532 nm was 
used to drive the experiment. The OPAs are of a bow- 
tie cavity design with a periodically poled KTP crystal 
providing the nonlinear interaction. In order to max- 
imise the cavity escape efficiency, faces of the nonlinear 
crystal were anti-reflection coated using the ion beam 
sputtering technique providing ultra-low loss for the 
intra-cavity circulating field. Also the cavity mirrors 
were custom made with reflectivities exceeding 99.95 %. 
The squeezed field is exiting the cavity through an out- 
put coupler, which was of 10 % transmissivity. The 
two squeezed beams, with almost perfectly equal per- 
formance with squeezing levels of approximately -6 dB 
and an anti-squeezing of 8.5 dB, are mixed on a 50:50 
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Figure 2: An experimental schematics used to generate and 
analyse a pair of entangled beams. OPA: optical parametric 
amplifier, BS: 50:50 beamsplitter, LO: local oscillator, HD: ho- 
modyne detection, (p: relative phase between optical fields. 



beamsplitter and locked in quadrature; this is keeping 
their relative phase shift (pent of re/ 2. This results in 
the generation of a pair of entangled beams, which are 
analysed using two balanced homodyne detectors (HD). 
These detectors had a quantum efficiency of 95 % and a 
dark noise clearance of 20 dB at the sideband detection 
frequency of 3 MHz. For the single photon reconstruc- 
tion we used very weak squeezing of 0.8 dB at a side- 
band frequency of 6 MHz. 

In order to characterise the Gaussian entanglement, 
we need to successively measure second order corre- 
lations and anti-correlations between the HD signals, 
when the HD systems are locked to the phase and 
amplitude quadratures, respectively. Measurements of 
both quadratures are accessible for the two entangled 
beams depending on the relative phase between the lo- 
cal oscillator (LO) beam and the analysed field. Hence 
we locked both HD systems to a half- or full-fringe of 
the interference pattern on the HD beamsplitters in or- 
der to measure the phase and amplitude quadratures, 
respectively. 

The homodyne detectors' AC and DC electronic sig- 
nals are recorded by a high speed data acquisition sys- 
tem (DAQ) and are used either to monitor the optical 
phase of individual beams in real time or for measure- 
ments of quantum correlations between the entangled 
fields. In particular, the phase lock between the two 
squeezed fields interfering on the entanglement beam 
splitter (BS) is crucial since it influences the quality of 
entanglement between the two BS output fields. The er- 
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ror signal for this lock was generated by a combination 
of HD detectors' DC signals using the DAQ. However, 
the error signals' DC offset fluctuations, having origin in 
parasitic interference, were setting the limit for the qual- 
ity of the entangled fields. We therefore corrected the er- 
ror signal for these fluctuations using the demodulated 
HD detectors' AC signals at a specific modulation fre- 
quency The corresponding arithmetics and the control 
loop itself were implemented in an FPGA. As a result a 
strong and stable error signal was produced allowing us 
to access states of high degree of entanglement. 

We evaluated several entanglement criteria based 
on the measured data. For the Duan inseparability 
criterion ll25il26l : 

^ = i ((^ (x« + x^)') + {Pa - P^)')) , (1) 

we found a value of I = 0.3 ± 0.02. We also evaluated 
the Reid and Drumtnond EPR criterion |[l9l |27| based on 
conditional variances: 

S = A^Klb^^b (2) 

with A^X^i^ denoting the conditional variance. Our data 
showed £ = 0.36 ± 0.03, which displays the lowest 
value reported so far to our knowledge. 

Single photon detection and evaluation scheme. The hy- 
brid experiments involving a single photon detector in 
one of the pair of modes take the homodyne detector 
data of the other mode into account only when the sin- 
gle photon detector fires. For the reconstruction of the 
quadrature histograms this means, that most of the ho- 
modyne data contributes with a weight of to the his- 
togram. Only when the single photon detector fires, the 
weight for the homodyne data is changed to 1. 

A pure two mode squeezed state can be written in the 
Fock basis as: 

-| 00 

= E(-tanhr)" |n), \n),, (3) 

with a and b for the two modes and r for the squeezing 
parameter. Hence, when n photons are found in mode a, 
the other mode b also has n photons. This is how most 
heralded single photon sources work. There is a click- 
detector in one mode. A click usually occurs, when there 
was at least one photon hitting the detector. So with ev- 
ery click in one mode there is at least one photon in the 
other mode. Now there is a tradeoff for the magnitude 
of r. On one hand the maximum probability of \ for one 
photon in each mode occurs for r = arcoshy^. How- 
ever, the probability for having two or more photons in 
this case is just as high, which leads to a highly mixed 
state. On the other hand for r ^ 0.1 the probability for 
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Figure 3: Quadrature probability distributions. The thin green 
curve represents the vacuum state. The thick blue curve shows 
the average of all quadrature measurements without condi- 
tioning. We carefully checked that this is independent of the 
projection angle. This is also true for the conditioned data, as 
shown by the thin lines in the background. The noise is due 
to statistics. The average of these curves is shown as thick red 
line. The dashed grey line indicates the result of an adjustment 
for the finite quantum efficiency. 

one photon drops to 1 % but two or more photons now 
occur one hundred times less likely. 

We extend the same idea with a continuous weight- 
ing based on data from the homodyne detector, which 
replaces the single photon detector in our case. 

(n) = i(x2 + p2-2) 

= i((Qi-i)),, (4) 

with 0^ denoting the average over all detection phases 
e and Qo = cos{e)X + sin(^)P. 

For the ordinary single mode tomography the prob- 
ability distribution for a given quadrature is recon- 
structed by taking statistically significant number of 
measurements cja of the same quadrature. The range of 
all possible measurement results is then divided into a 
certain number of discrete bins. Each time a measure- 
ment result falls into a particular bin, its value is in- 
creased by a fixed increment. This means that every sin- 
gle measurement represents the ensemble equally well. 
In our case the homodyne measurement result cji^ of the 
second mode b of the entangled pair provides some in- 
formation about the first mode a. When building the 
histogram for mode a instead of using a fixed incre- 
ment, the bin corresponding to the measurement qa is 
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increased by the increment df which depends on qi,: 

di = ql-\. (5) 

Like for ordinary tomography each of the detection 
phases is measured repeatedly to provide statistical sig- 
nificance. However, during the repeated measurements 
of mode a the detection quadrature of mode h is varied 
such that all phases occur uniformly According to Eqn. 
11 this means that the reconstructed histograms fhi^a) 
give the probability distributions of convoluted by 



(6) 



where is the mean value conditioned to Qa = c]a- 

We now note that for small squeezing parameters r, the 
probability for having 2 or more photons is small so 
we can approximate (n) ^ p{l) and hence {'nh)\qa — 
p{ljj\qa)- Injecting this result into ^ and using Bayes 
law to 'invert' conditional probabilities, we find 



(7) 



Within a renormalisation factor p(l^), we find the de- 
sired probability distribution corresponding ideally to a 
single photon in a heralded by a 'click' in b. 
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Figure 4: Reconstructed photon number distributions. The 
blue bars are extracted from the measurement data directly. 
The red bars are the result of the single photon conditioning. 
The single photon contribution exceeds 67%, which we would 
like to consider as a benchmark figure for the overall quality 
of the experimental apparatus. Note that there is no correction 
for detection efficiency or dark-noise for these data. The grey 
bar indicates the result of an adjustment for the finite quantum 
efficiency. 

The complete characterisation of a general two mode 
state by means of quantum state tomography requires 



the measurement of a large number of projection an- 
gles and their combinations, which is experimentally 
challenging. In the particular scenario considered here, 
the requirements are far less complex. The individual 
modes of an entangled two mode squeezed state are 
known to be circular symmetric in quadrature phase 
space. However, the correlation between the two de- 
pends on the relative detection phase. The condition- 
ing method we use for the single photon projection av- 
erages over all detection phases for one of the modes. 
This implies, that also the dependance of the mutual de- 
tection phase is cancelled out. For the reconstruction we 
ensured, that there actually is no phase dependence in 
our data. This allows us to base our tomographic recon- 
struction on a single averaged quadrature distribution, 
which we assume to be the same for all quadratures. 
Figure [3] shows that this assumption is valid. 

On the experimental side we ensured that for the de- 
tection of the conditioning mode b all quadrature pro- 
jections are represented equally by rapidly scanning the 
optical phase of the local oscillator with a range adjusted 
to cover an interval of In. At the same time we slowly 
varied the detection phase of the detection of mode a 
over a range much wider than In. The data collection 
time for each histogram was chosen to be short com- 
pared to the variation of mode a and also as an integer 
multiple of the scanning period of mode b. That means 
that each histogram represents a fixed quadrature pro- 
jection for mode a with a uniform contribution of all pro- 
jection phases of mode b. From this procedure we see, 
that the shape of the reconstructed probability distribu- 
tion indeed does not depend on the detection phase of 
mode a. Thus, we were able do the state tomography as- 
suming a circular symmetric state without exactly keep- 
ing track of the detection phase. 

Based on this probability distribution we used the in- 
verse radon transform technique for the reconstruction 
of the Wigner function, see Figure [l] Furthermore we 
used the pattern function method of |i28J in order to re- 
construct the Fock base density matrix. The diagonal 
elements are shown in Figure [4] Clearly, the n = 1 
component is the most dominant one. The residual 
vacuum component arises from experimental imperfec- 
tions, such as preparation or detection efficiencies. For a 
perfect experiment the vacuum component would van- 
ish. We used the capability of the reconstruction method 
to adjust for the finite efficiency for the grey components 
of Figure |4] and |3]. In this sense we consider the residual 
vacuum component as an overall figure of merit for the 
quality of the experiment. 

Conclusion. The advances in technology and meth- 
ods enable the implementation of better and more elab- 
orate quantum optics experiments. In our experiment 
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this gives rise to a remarkable progress (see Figure 9 
in 1 20 1) for the violation of the Reid and Drummond 
EPR criterion as well as the Duan inseparability crite- 
rion. Furthermore, sophisticated data acquisition and 
processing techniques allowed us to apply the proposal 
of [21Jand reconstruct a high purity single photon state 
from our measurement data. This complementary ap- 
proach to single photon detection could be useful for 
the implementation of more complex protocols which 
might be prohibited by the limitations of single photon 
detection. 
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